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Abstract 

In a previous publication, we have constructed a set of non-linear evolution equations for dipole 
scattering amplitudes in QCD at high energy, which extends the Balitsky-JIMWLK hierarchy by 
including the effects of fluctuations in the gluon number in the target wavefunction. In doing so, 
we have relied on the color dipole picture, valid in the limit where the number of colors is large, 
and we have made some further approximations on the relation between scattering amplitudes 
and dipole densities, which amount to neglecting the non-locality of the two-gluon exchanges. 
In this Letter, we relax the latter approximations, and thus restore the correct structure of the 
'triple-pomeron vertex' which describes the splitting of one BFKL pomeron into two within the 
terms responsible for fluctuations. The ensuing triple-pomeron vertex coincides with the one 
previously derived by Braun and Vacca within perturbative QCD. The evolution equations can 
be recast in a Langevin form, but with a multivariable noise term with off-diagonal correlations. 
Our equations are shown to be equivalent with the modified version of the JIMWLK equation 
recently proposed by Mueller, Shoshi, and Wong. 
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1 Introduction 



In a recent publication [1], we have shown that the Balitsky-JIMWLK (Jalilian- 
Marian-Iancu-McLerran-Weigert-Leonidov-Kovner) hierarchy of equations [2-5] previ- 
ously proposed to describe the non-linear evolution of scattering amplitudes in QCD at 
high energy fails to include the effects of fluctuations in the gluon number in the target 
wavefunction, that were recently demonstrated to play an important role in the approach 
towards saturation and the unitarity limit [6, 7]. Still in Ref. [1], we have then constructed 
a generalization of these equations which includes the relevant fluctuations in the limit 
where the number of colors is large (N c ^> 1) and in the kinematical region where these 
effects are expected to be important [6,7], namely, in the low density regime at high 
transverse momenta (well above the target saturation momentum [8,9]). To that aim, we 
have relied on Mueller's color dipole picture [10-12], which is indeed known to provide 
a complete description of the target wavefunction (including the relevant fluctuations) 
in the approximation of interest. Besides, in order to simplify some of the technical ma- 
nipulations in Ref. [1], we have performed some further approximations in the relation 
between scattering amplitudes and dipole densities, which are tantamount to neglecting 
the non-locality of the two-gluon exchanges in dipole-dipole scattering. 

Although reasonable within the most interesting kinematical region — namely, within 
the intermediate regime at high momenta which is characterized by the BFKL 'anoma- 
lous dimension' [8, 13-17] — , these additional approximations have ineluctably entailed 
a loss of accuracy, and also some conceptual shortcomings, for our final results in Ref. 
[1]. Specifically, the normalization of the fluctuation terms in our final equations there is 
afflicted with an unknown 'fudge' factor, and the non-locality structure of those terms is 
oversimplified. Because of that, these equations cannot be used, e.g., to study the tran- 
sition towards the high-Q 2 evolution [18], which is expected to take place at sufficiently 
large transverse momenta. For instance, one would like to understand the interplay be- 
tween the fluctuation-dominated dynamics predicted by the analysis in Refs. [1, 7] (which 
should completely supersede the BFKL dynamics at very high energies) and the DGLAP- 
like dynamics [18] which should emerge at higher transverse momenta. But in order to 
do so, one would need the correct non-local structure of all the terms in the equations, 
including the fluctuation terms. Besides, the too drastic simplifications performed on the 
structure of the fluctuation terms leads to numerical difficulties which prevent a direct 
resolution of the hierarchy of equations proposed in Ref. [1] through iterations [19]. These 
difficulties are eliminated by an additional coarse-graining in impact parameter space, 
which leads to a mathematically well-defined Langevin equation [1] - the QCD analog 
of the stochastic Fisher-Kolmogorov-Petrovsky-Piscounov (sFKPP) equation [20]. But 
this final equation has the further drawback to ignore the impact-parameter dependence 
of the scattering amplitudes. Last but not least, we would like to be able to compare 
the newly proposed terms in Ref. [1] with related results in perturbative QCD [21-26]. 
As explained in [1], these terms describe pomeron splittings in the target wavefunction. 
The relevant triple-pomeron vertex has been computed in perturbative QCD in the re- 
cent years [21,23,25,26], and we would like to identify this as a building block in our 
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fluctuation terms. 



For all these reasons it would be important to relax the additional approximations 
performed in Ref. [1], and thus obtain the complete equations valid to leading order at 
large N c . This is the task that we shall undertake in this Letter. Specifically, by using the 
complete relations between scattering amplitudes and dipole densities valid in the dilute 
regime, we shall be able to restore the full non-locality of the fluctuation terms in the 
evolution equations for the amplitudes, and thus correct all the deficiencies alluded to 
before. In particular, in our new equations we shall recognize, as expected, the large-iV c 
version of the triple-pomeron vertex describing pomeron splitting, in the form given by 
Braun and Vacca [23] (see also [22]). Since we knew already that the triple-pomeron 
vertex responsible for pomeron mergings is also included [25,27], through the non-linear 
terms in our equations (which are the same as the corresponding terms in the Balitsky- 
JIMWLK hierarchy [2-5]), we conclude that the new equations generate correctly, through 
iterations, the pomeron loops expected in perturbative QCD at large N c [11, 12]. 

In view of numerical studies, it is often convenient to dispose of a stochastic formula- 
tion of a hierarchy of equations. We shall show that our new equations can be equivalently 
reformulated as a Langevin equation, which generalizes the one proposed in Ref. [1] by 
including the impact-parameter dependence and the full non-locality of the coefficient 
of the noise term. The counterpart of this is the fact that the new noise term depends 
upon three variables, and has non-diagonal correlations in two of them, which rises some 
doubts about his usefulness for practical calculations. 

In this paper, we shall follow the same general strategy as in Ref. [1], namely, we shall 
view the high-energy evolution as the evolution of the target wavefunction, for which we 
shall use the dipole picture [10] in the dilute regime at high momenta and the JIMWLK 
evolution [3-5] in the high density regime near, or at, saturation. Very recently, Mueller, 
Shoshi, and Wong [28] have proposed a generalization of the JIMWLK equation which 
includes the effects of pomeron splitting in the dilute regime, at the same accuracy level as 
the dipole picture. We shall explicitly check that our new equations are indeed equivalent 
with the evolution described in Ref. [28]. In particular, the need for a non-diagonal noise 
term in order to achieve a stochastic formulation has been originally recognized in [28]. 

The paper is structured as follows: In Sect. 2, we shall construct the fluctuation term 
for the equation describing the evolution of the scattering amplitude for a pair of dipoles 
(which is the lowest equation in the hierarchy where this term enters). Then, in Sect. 
3, we shall generalize this construction to the equations of higher level in the hierarchy 
(pertinent to the scattering of k dipoles with k > 3), and we shall introduce the alternative 
formulation of the hierarchy as a Langevin equation. In Sect. 4 we shall demonstrate the 
equivalence with the extended JIMWLK formalism of Ref. [28]. Finally, in Sect. 5 we shall 
analyze the structure of the triple-pomeron vertex encoded in our equation, and recognize 
this as the vertex originally computed by Braun and Vacca [23]. 



3 



2 The scattering amplitude for a pair of dipoles 



In this section we will derive the evolution equation for the amplitude describing the 
scattering of a dipole pair off an arbitrary hadronic target in the large-iV c limit. We shall 
view the scattering in a frame in which most of the total rapidity Y, with a s Y > 1, 
is carried by the target (whose wavefunction is therefore highly evolved), whereas the 
projectile is simply a set of elementary dipoles. Also, we shall discuss explicitly only the 
contributions to the evolution arising from fluctuations in the gluon number in the target 
wavefunction. From the discussion in Refs. [1,6,7], we expect such fluctuations to play 
an important role only in the tail of the distribution at large transverse momenta, where 
the gluon occupation factor is low (of order one), and the dipole picture is reliable (for 
large N c ). To study this low density region, we shall assume that the external dipoles 
composing the projectile have small transverse sizes (much smaller than the inverse of the 
target saturation momentum at their respective impact parameters). 

Consider first the scattering of a single external dipole, with transverse coordinates 
x and y (for the quark and antiquark, respectively). In the present approximations, this 
dipole undergoes single scattering via two gluon exchange with any of the dipoles internal 
to the target. Thus, the corresponding scattering amplitude reads: 

(T(x,y)) Y = a 2 J d 2 ud 2 vA (x,y\u,v)n Y (u,v) (2.1) 

where, up to a factor a 2 , Ao is the amplitude for dipole-dipole scattering in the two-gluon 
exchange approximation and for large N c : 



A (x,y\u,v) = i 



' (x - v) 2 (y - u) 2 
(x — u) 2 (y — v) 2 



; 2 . 2 ) 



The factor of a 2 has been extracted from the amplitude but included explicitly in Eq. (2.1), 
in order to facilitate the perturbative counting for the subsequent expressions. Further- 
more, n Y (u,v) measures the average number of dipoles with coordinates (u,v) in the 
wavefunction of the hadron as produced after a BFKL evolution over the rapidity Y. 

At this point, and for later convenience, it is useful to invert Eq. (2.1) and express 
the dipole density in terms of the scattering amplitude. Applying V 2 x V 2 y on Eq. (2.1) and 
using the property V^lna? 2 = Att5^ 2 \x), we obtain 

n Y (x, y) + n Y (y, x) = ^V 2 x V 2 y (T(x, y)) Y , (2.3) 

with the assumption x ^ y. 

Now let us consider the scattering of a pair of dipoles (xi,yi) and (012,2/2) °^ ^ ne 
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target. Clearly the extension of Eq. (2.1) reads 



(T {2 \x 1 ,y 1 ;x 2 ,y 2 )) Y = a t J A(^i, «i) Ao(x 2 , y 2 \u 2: v 2 ) n { y\u 1: v ± ; u 2 , v 2 ), 

Ui,Vi 

(2.4) 

where rty\ui, v±; u 2 , v 2 ) is the average number density of pairs of dipoles in the target. 
As in Ref. [1], we include in n Y only the pairs made of two different dipoles, that is, 
we do not allow the two external dipoles to scatter off a same dipole in the target. This 
should be a good approximation so long as the external dipoles have different sizes and/or 
impact parameters. In particular, if at Y = the target consists in a single dipole (this 
will be our initial condition), then = 0, and therefore (T^} = as well. 

The evolution equation for (T^) Y can be obtained from the corresponding equation 
for , which has been derived in Ref. [1]. The inhomogeneous part of the latter, i.e. the 
piece which describes low-density fluctuations, reads 



= ^M(u l ,v 2 ,u 2 )n Y (u 1 ,v 2 )5 (2 \u 2 -v l ) + {l <-> 2}, (2.5) 



fluct 27r 



dY 

with the dipole kernel [10] 

M(x,y,z) = - {X ~ y)2 (2.6) 

(x - z) 2 (y - z) 2 

(The terms not shown here correspond to the usual BFKL evolution of the dipole pair 
density; see Eq. (5.15) in Ref. [1].) Using Eqs. (2.4) and (2.5) to write the evolution of 
(T( 2 )) y in terms of the dipole density n Y and then employing Eq. (2.3) in order to express 
n Y in terms of the scattering amplitude (T) Y , it is straightforward to show that 



dY 



OL \ 2 OL f 

2tt) 2n J M ( u > v > z ) A o( x ±>y±\ u > z ) A °( X2 >y 2 \ z > v ) 

xVX (T(u,v)) Y . (2.7) 



fluct u,v,z 



The complete evolution equation for (T^)y can be now obtained by restoring the 
more standard terms associated with the BFKL evolution and with unitarity corrections 
at high energy. The former are terms linear in (T^)y, while the latter are non-linear 
terms proportional to the scattering amplitude (T^) Y for a system of three external 
dipoles, and are the same as in the corresponding Balitsky-JIMWLK equation [2, 4, 5] 
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(here, for large N c ). Thus our final evolution equation for (T^)y reads 



d(T( 2 \x 1 ,y 1 ;x 2 ,y 2 )^ 



Y 



<)) = £ / { i MxiVlZ ^ ( T(2) ( Xl ' 2/15 X2 ' ^ 

z 

-A<(aJi, (T( 3 )(a;i,z;«,yi;a! 2 ,i/ 2 )) y ] + [1 <-> 2]} 

d(TW(x 1 ,y 1 ;x 2 ,y 2 ) 



+ ^ ^ 



(2.J 

fluct 



<9F 

with the last term given by Eq. (2.7) and where we have introduced the shorthand notation 

M xyz <g> f(x, y) = M{x, y, z)[-f(x, y) + /(*, z) + f(z, y)\. (2.9) 

Under the same approximations, the scattering amplitude for a single dipole obeys the 
following equation, which is formally the same as the first equation in Balitsky hierarchy: 

d{n * Y y)) Y =|/ [M xyz ® (T(x, y) Y - M(x, y, z) (T^(x, z; z, y)) y }. (2.10) 



z 



Note that the evolutions of (T) Y and (T^)y are coupled by the fluctuation term in 
Eq. (2.8), in agreement with [1]. Moreover, (T^ 2 ^)y is also coupled to (T^)y, by the 
unitarity corrections, so the equations above do not close by themselves, but rather are a 
part of an extended (actually, infinite) hierarchy. The higher equations in this hierarchy 
will be written down in the next section. 

Before closing this section, we would like to comment on the simplified version of 
Eq. (2.7) which was proposed in [1]. There, an approximation was done at the level of 
Eq. (2.2), in which the non-locality of the 2-gluon exchange interaction was neglected. 
That is, the elementary dipole-dipole scattering amplitude was effectively assumed to be 

a 2 s A (x, y\u, v) ~ ca 2 s (x - y) 4 S {2) (x - u) 5 {2 \y - v), (2.11) 

where c is an unknown "fudge" factor. Then, by following the same steps as above, it was 
shown that the right hand side of Eq. (2.7) is 



After a coarse-graining in impact parameter space and a saddle point approximation, 
a Langevin equation, equivalent to the stochastic FKPP equation, was finally obtained. 
However, these equations suffer from the shortcomings mentioned in the Introduction. 
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3 The Langevin Equation 



In order to have a complete hierarchy of equations, we need to generalize Eq. (2.7) 
to the case where an arbitrary number of dipoles, say k, scatters of the target. Given 
the K-th density of dipoles ny(ui, V\\ u K , v K ) = ny({u, v}) in the wavefunction of the 
hadron, the amplitude for K-dipoles scattering reads 

T iK \{x, y})) y = a 2 s K J Ao(x u i/i|ui, Vl ) ... A (x K , y K \u K , v K ) n<y\{u, v}). (3.1) 

Ui,Vi 

The corresponding evolution equation is obtained by the straightforward generalization 
of the procedure that we followed in the previous section. The rate of change, due to 
fluctuations, of the K-th density of dipoles of the hadronic wavefunction depends only in 
the (k — l)-th density, and therefore the same will be true for the scattering amplitude 
(T^ K \{x, y}))y, since Eq. (3.1) relates quantities of the "same order" in k. The ensuing 
equation reads 

d(T^{{x,y})\ fa s \ 2 a s f 

dY ~ = \2ttJ 27r J M ( u ' v ' z "> A o(x 1 ,y 1 \u,z)A (x2,y 2 \z,v) 

fluct u,v,z 

xVX (T^- 1 \u,v;x 3 ,y 3 ;...;x K ,y K )) Y + perm., 

(3.2) 

where "perm." stands for all possible permutations in the arguments of (T( K_1 ))y, so that 
the total number of terms in the right hand side is k(k — l)/2. 

It is now straightforward to add to the r.h.s. of Eq. (3.2) the (large- iV c version of the) 
standard Balitsky-JIMWLK terms, and thus obtain the complete hierarchy for (T^) Y 
within the present approximation. Since in practice (e.g., for numerical studies) it seems 
difficult to deal with a hierarchy as a whole, it is convenient to notice that one can 
equivalently reformulate this hierarchy as a single stochastic equation. Let us consider the 
following Langevin equation 



dT Y (x,y) 
dY 



fluct 



= 271/27 J MX ' ylu > Z) u - z) 2 TAUj V) V[U ' Vj Z ' Y) ' 



(u-zf 

(3.3) 



where the noise satisfies 



(u(u u Vl , z u Y)u(u 2 , v 2 , Z2, Y')) = 5^ 2 \ Ul - v 2 )5^( Vl - u 2 )S^( Zl - z 2 )5(Y - Y') 

(3.4) 

with all the other correlators being zero. Note that the correlator (3.4) is non-diagonal in 
the first two arguments (it, v) of the noise. A similar property has been recently noticed 
in Ref. [28], in a physical context which is very close to ours (see the discussion in Sect. 4 
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below), but we are not aware about previous studies of such a non-diagonal noise in other 
(physical or mathematical) contexts. The stochastic problem above must be understood 
with the Ito prescription for rapidity discretization (see, e.g., [1,29] for details). Then 
one can check that by solving this equation one generates the same correlations as from 
the original hierarchy in Eq. (3.2). Indeed, consider the change in the K-dipole function 
T\...T K under an increase AY in rapidity (in Tj, the index % serves as a shorthand for 
the coordinates of the i-th dipole). Eqs. (3.3) and (3.4) imply that AT ~ V AY, so in 
evaluating A(T 1 ...T /t ), it is enough to keep terms of second order in the AT's : 

A(Tl..T k ) = (7\ + ATx)...(T K + AT K ) - T X ...T K - AT, AT 2 T 3 ... T K + perm. (3.5) 

Indeed, the linear terms will vanish after averaging over the noise, whereas the cubic 
terms, or higher, will be of higher order in AY, and thus will vanish in the continuum 
limit Ay — > 0. The transverse coordinates dependence of the noise correlation in Eq. (3.4) 
is such that, when we take the expectation value in Eq. (3.5), the emerging hierarchy is 
exactly the same as the one in Eq. (3.2), provided we make the identification (Ti...T K ) <-> 
(T^({x,y})) Y . 

Eq. (3.3) refers solely to the fluctuations terms in the hierarchy, but the standard 
Balitsky-JIMWLK terms can be easily added [1]. We finally obtain 



dT Y (x,y) _a ± 
dY ~ 2tt 

+ 



J [M xyz <g> T Y (x, y) - M(x, y, z) T Y (x, z) T Y (z, y 
dT Y (x,y) 



dY 



(3.6) 

fluct 



with the last term being given by Eq. (3.3). This equation together with the noise term 
(3.4) is the Langevin problem equivalent to the non-linear evolution equations with 
pomeron loops. But, as mentioned above, the non-diagonal nature of the noise corre- 
lator (3.4) is a rather unusual mathematical property, and presently we do not know 
whether such a noise can be useful for numerical studies, or not. 



4 The Color Glass Condensate Approach 



An alternative approach to high-energy QCD is the Color Glass Condensate (CGC) 
formalism in which one studies the random color fields created by 'color sources' (gluons 
or dipoles) within the target, and the evolution of the respective correlation functions 
with increasing Y (see [30] for review papers and more references). In a specific gauge, 
the color field in the target has only one component, the light-cone component A£, which 
is traditionally denoted as A+(x) = a a (x~, x). (This is a function of the longitudinal coor- 
dinate x~ together with the transverse coordinate x.) Let WV[q] represent the functional 
probability distribution for finding a particular configuration of this field at rapidity Y. 
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Then the expectation value of an operator O[o\ is given by 



(0[a]) Y = J 



DaW Y [ot]0[a\. 



(4.1) 



For example the operator for the scattering between a single dipole (x, y) and the CGC 
is computed in the eikonal approximation as 



where is a Wilson line describing the scattering between the quark and the color field 
in the CGC (below, V stands for a path ordered product), 



while V y refers similarly to the antiquark. Furthermore, the corresponding operator for a 
dipole pair is T^(x 1 ,y 1 ;x 2 ,y 2 ) = T(x 1: y 1 )T(x 1 ,y 1 ), with T given by Eq. (4.2). 

Under an increase dY in rapidity, the correlation functions of the color field change 
because of the induced emission of a small-a; gluon. These changes can be accommodated 
in a functional evolution equation for the probability distribution, which can be written 
in Hamiltonian form (at least, formally) : 



Then the evolution equation for (0[ai])y can be deduced from Eqs. (4.1) and (4.4). 

The original version of the Hamiltonian which enters Eq. (4.4), known as the JIMWLK 
Hamiltonian (-^jimwlk) [3-5], leads to Balitsky equations [2] when applied to scattering 
operators like (4.2). -£/jimwlk is a second-order functional differential operator, with a 
kernel which is non-linear in the field a to all orders. The quadratic part of this kernel is 
responsible for the linear, BFKL, evolution, while the terms of higher order in a generate 
non-linear terms, so like the terms expressing unitarity corrections in Balitsky equations. 
As we have noticed in [1], the evolution generated by -£/jimwlk does not include gluon 
splitting, and therefore cannot describe the correlations associated with fluctuations in the 
gluon number, which are however important in the dilute regime. Very recently, Mueller, 
Shoshi and Wong [28] have proposed an additional term Ha nct , which is inspired by the 
probabilistic interpretation of the evolution in the dipole model [10,31], and which takes 
into account the low density fluctuations at the same level of accuracy as the dipole 
picture. Thus, the total Hamiltonian used in Ref. [28] can be written as: 




(4.2) 




(4.3) 



_d_ 

dY 



W Y [a] =HW Y [a]. 



(4.4) 



H — H JIMWLK + -ftfluct, 



(4.5) 
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where -Hjimwlk can be found in Refs. [4,5] (see also [32] for a simple expression), and 



Ha uct = ^ J M(u,v, z)g(u l \u, z)g(v 1 \u, z)g(u 2 \z,v)g(v 2 \z,v) 

x f 7 n f 7 n f U n7 J , V*Vfc»c*». (4.6) 
5a a (wi) 5a a («i) fe fc (w 2 ) 5a fe (v2) 

Here, the integration goes over all the transverse coordinates u, v, z, Ui, v±, u 2 , v 2 , and 
a a (x) = / rfx~o; a (a; _ , x). The function Q(ui\u, z) is, up to a factor gt a , the classical field 
created by the elementary dipole (u,z), and reads 

g(u 1 \u,z) = ^-\n ( ; Ul ~ z) * (4.7) 

47T (Wi - w) 2 



In what follows we would like to show that the Hamiltonian (4.6) leads precisely to our 
previous hierarchy (3.2) for the dipole scattering amplitudes. We shall demonstrate that 
only for the case of two external dipoles (k — 2), since the corresponding proof for arbitrary 
k follows trivially. Also we shall focus on the contribution induced by the fluctuation piece 
.f/fiuct of the Hamiltonian, since the corresponding contribution of //jimwlk is well-known 
[4, 5, 29, 32] to be the respective Balitsky equation. Thus, we start with (cf. Eqs. (4.1) and 
(4.4)): 



d(TW( Xl , yi ;x 2 ,y 2 
dY 



fluct 



= J ' DaT^(x 1 ,y 1 ;x 2 ,y 2 )H &uct W Y [a]. (4.8) 



In the weak field limit the scattering amplitude is determined from the small coupling 
expansion of the corresponding Wilson lines operator and it reads 

T&>( Xl , y i; x 2 , y 2 ) = [a d ( Xl ) - a d ( yi )] 2 [a e (x 2 ) - a%y 2 )] 2 . (4.9) 

By partial integrations within the r.h.s. of Eq. (4.8), one can bring the functional differen- 
tial operator (5 /5a) 4 (cf. Eq. (4.6)) to act on T^ 2 \ This operation annihilates the a fields 
manifest in Eq. (4.9) and gives rise to a sum of products of two-dimensional delta func- 
tions. For consistency with the derivation of Eq. (4.6) [28], we should restrict ourselves to 
the terms of leading order in iV c , which arise when a pair of derivatives which are traced 
over color acts on a pair of fields which carry both the same color indices. Then we can 
perform the integrations over u 1 ,Vi,u 2 ,v 2 , to finally arrive at 



d(T^(x 1 ,y 1 ;x 2 ,y 2 )) Y «3 n 



dY 



fluct c u>Vt 



J M(u, v, z) A (x u j/xlu, z) A (x 2 , y 2 \z, v) 
xVX (-2a c (u)a c (v)) Y . (4.10) 
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The last factor in the above equation may be rewritten as 

V^V 2 (-2a%u)a%v)) Y = V^V 2 ([a c (u) - a%v)]% = ^ y^V 2 (T(u, v)) Y , 

(4.11) 

where we have also used the weak-field expansion of the scattering amplitude (4.2) for a 
single dipole. With this identification, Eq. (4.10) leads indeed to Eq. (2.7), which completes 
our proof. 



5 The triple pomeron vertex 

The above derivation of our new equations within the CGC approach sheds more light 
on their physical interpretation and confirms the original discussion in Ref . [1] : when the 
evolution is included in the target wavefunction, the fluctuation terms in the equations 
for the amplitudes (e.g., the term linear in (T^ k ~^)y in the r.h.s. of Eq. (3.2) for (T^)y) 
correspond to gluon splittings, by which we mean more precisely the transition vertex from 
two to four gluons described by the new piece H^ nct , Eq. (4.6), of the CGC Hamiltonian. 
Similarly, the non-linear terms responsible for unitarization (e.g., the term (T^)y in the 
r.h.s. of Eq. (2.8) for (T^)y), correspond to gluon mergings, as described by the terms 
of higher order in a in the original JIMWLK Hamiltonian. By 'gluon mergings' we mean 
all the vertices leading to a reduction in the number of gluon fields. The simplest among 
them is the transition from four to two gluons, and is described by that term in -ffjiMWLK 
which is dual to Eq. (4.6), in the sense of containing four powers of a and two powers 
of the derivative {5 /5a). In general, one would expect this duality to extend to higher 
vertices; that is, H^ uct is expected to include terms with higher-order derivatives, which 
should be 'dual' to the terms with higher powers of a in -^jimwlk- Since such terms do not 
appear to be generated by the dipole picture which has been used to study fluctuations 
so far [1, 6, 7, 28], they are presumably of higher order in 1/N C . 

Furthermore, within the framework of the large- N c approximation, the gluon number 
changing transitions alluded to above should be tantamount to pomeron splittings, or 
mergings, i.e., processes through which one BFKL pomeron is connected to two. For 
the non-linear term in the Balitsky-Kovchegov equation [2, 33] (and thus, by extension, 
for the whole Balitsky hierarchy at large N c ), this has been checked already in Refs. 
[25, 27]. Namely, it has been shown there that the general 2-to-4 gluon vertex [21] reduces 
essentially to the dipole kernel, Eq. (2.6), when restricted to the large-iV c limit and to 
the physical conditions appropriate for the merging of two pomerons 2 . This is indeed 
consistent with the structure of the non-linear term in Eq. (3.6). 

2 Namely, the four gluons on one side of the vertex are pairwise connected in color singlet 
states, which moreover are multiplied by functions which vanish at equal coordinates, so like 
T Y (x,z)T Y {z,y) in the r.h.s. of Eq. (3.6). 
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In what follows, we shall present a similar check for the vertex that enters the fluc- 
tuation term in Eq. (2.7) (and thus the noise term in Eq. (3.6)). To that aim, we shall 
compare our result for this vertex with that in Ref. [23], where a linear evolution equation 
has been constructed which includes the triple-pomeron vertex for one pomeron splitting 
into two. (The equation written in Ref. [23] refers to the diffractive amplitude for four 
gluons, which is quite similar to the scattering amplitude for two dipoles in our approach.) 

From our equations, the amputated version of the fluctuation vertex can be imme- 
diately extracted, as the function which is convoluted with (T(u,v)) Y in the r.h.s. of 
Eq. (2.7). But since in Ref. [23] one rather finds the non-amputated version of the triple- 
pomeron vertex (this is referred there as the 'triple pomeron interaction', or TPI), we 
find it convenient to compute the similar quantity in our approach as well, by solving the 
linearized version of Eq. (2.8). Assuming the initial condition (T^) = 0, the solution to 
this inhomogeneous equation can be written as 

Y 



(T {2) (x u y 1 ;x 2 ,y 2 )) Y = J dy J G Y - y (x 1: yi\u u v ± ) G Y _ y (x 2 , y 2 \u 2 , v 2 ) 

d(T (2 \u l ,v 1 ;u 2 ,v 2 



Ui,Vi 

d(T^(u 1 ,v 1 ;u 2 ,v 2 )) 

(5.1) 

d V fluct 

where the Green's function G Y satisfies the BFKL equation (which is homogeneous) 

W ~^I Mxyz ® G Y (x : y\u,v) = : (5.2) 

z 

along with the initial condition 

G Y {x, y\u, v) = 5 {2 \x - u)5 i2 \y - v). (5.3) 

Instead of directly solving Eqs. (5.2) and (5.3), which is a simple task in any case, it is more 
convenient for our purposes to observe that one can combine each of the Green's function 
in Eq. (5.1) with one of the elementary amplitudes Aq contained in the source term 
d{T<®) y /dy, cf. Eq. (2.7). By doing this, we can reconstruct the dipole-dipole scattering 
amplitude in the BFKL approximation (the particular case of the amplitude (T(x, y)) Y in 
which the target is itself an elementary dipole). Indeed, since the dipole-dipole scattering 
amplitude in the linear regime satisfies the BFKL equation, one can write 

A Y - V {x 1 ,y 1 \u,z)= j G Y - y (x 1 ,y 1 \u 1 ,v 1 )A (u 1: y 1 \u,z) : (5.4) 

U1,V1 

along with an analogous equation for the u 2 ,v 2 integration. Therefore, we can finally 
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express the scattering amplitude for the dipole pair in the form 

2 Y 

(T ( - 2 \x 1 ,y 1 ;x 2 ,y 2 )) Y = (^j ^ j dy J M(u,v, z) Ay-y^y^u, z) 

u,v,z 

x A Y -v(x 2 , y 2 \z, v) V 2 u Vl (T(u, v)) y , (5.5) 

which is indeed recognized as the 'triple pomeron interaction' in [23] (see Eq. (69) there). 
This concludes our present analysis. 

Only a few days ago, when this work was already completed, a preprint appeared 
by E. Levin and M. Lublinsky [34] in which a hierarchy of equations including pomeron 
loops is proposed. Although the general structure of this hierarchy bears some similarity 
to ours, its derivation and detailed form remain different. 
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